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O ' Abstract 

Ph! We compute next-to-leading order QCD corrections to the CP asymmetry = Im(r 12 /M 12 ) 
q h " in flavour-specific B^s decays such as Bd — > XiT7 e or B s — > D~ir + . The corrections reduce the 
uncertainties associated with the choice of the renormalization scheme for the quark masses 
significantly. In the Standard Model we predict af s = —(5.0 ± 1.1) x 10~ 4 . As a by-product 
we also obtain the width difference in the B^ system at next-to-leading order in QCD. 
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1 Preliminaries 

Bd and B s mesons mix with their antiparticles. The time evolution of the B q — B q system 
(with q = d or s) is characterized by two hermitian 2x2 matrices, the mass matrix M q and 
the decay matrix T q . The oscillations between the flavour eigenstates B q and B q involve the 
three physical quantities |Mf 2 |, |rf 2 | and 4> q = arg(— Ml 2 /T\ 2 ) (see e.g. [T]). They are related 
to the mass and width differences of the B q system as 

AM, = 2\M? 2 \, AT q = r q L -r q H = 2|r? 2 |cos0 (? , (1) 
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Figure 1: Leading order contribution to Ti2 (left) and a sample NLO diagram (right). The 
crosses denote effective AB = 1 operators triggering the b decay. The full set of NLO diagrams 
can be found in [Zj. 



where T q L and T q H denote the widths of the lighter and heavier mass eigenstate, respectively. 
Here and in the following we neglect tiny corrections of order \T\ 2 /M- 



q i2 

12 1 • 



The CP-violating phase (J) q can be measured through the CP asymmetry a q s in flavour- 
specific B q — > / decays, which means that the decays B q — > / and B q —> f are forbidden j2]: 

-« - r «^ i >-«- r < B '< i )- 7 >-. m ^^ tm ,,. (2) 



UB,(t) - /) + r(B,(«) - /) Mh AM, 



'I 



Here B q (t) and B q (t) denote mesons which are tagged as a B q and B q at time t = 0, re- 
spectively. An additional requirement in Eq. (J2J) is the absence of direct CP violation in 
B q — > /, which is equivalent to \ (f\B q ) \ = \ (f\B q )\. For example, af s can be obtained through 
B s — > .D~7r + . The standard way to access a q s uses _B g — >• X£~Vi decays, which justifies the 
name semileptonic CP asymmetry for a q s . The measurement of a q s does not require tagging 
(see e.g. [3]). A further method to access a q s uses the fully inclusive, tagged B decay asymmetry 
discussed in jl]. 

a|, is small because of two suppression factors: First IT^/M^I = 0(m^/M^) suppresses a q is 
to the percent level. Second there is a GIM suppression factor w? c jm^ reducing a q is by another 
order of magnitude. This GIM suppression is lifted if new physics contributes to argMi2- 
Therefore aj s is very sensitive to new CP phases Up to now, the Standard Model (SM) 

prediction for a q s was only known in the leading-logarithmic approximation. The unknown 
next-to-leading order (NLO) QCD corrections were identified as the largest theoretical uncer- 
tainty in a q s [S]. While NLO corrections were calculated long ago for M q 2 [Hj, only certain 
portions of the QCD corrections to T\ 2 (relevant to Ar s ) were known so far In Sect. El we 
compute the missing pieces of the latter. Predictions for a\ s and Ar d can be found in Sect. El 



2 at next-to-leading order in QCD 

In this section we specify the discussion to the case q = d and omit the index q. The gene- 
ralization of our results to is straightforward. Ti2 is an inclusive quantity stemming from 
decays into final states common to B and B. It can be computed with the help of the heavy 
quark expansion (HQE) jHj from diagrams like those in Figure ^ The HQE is a simultaneous 
expansion in AQ CD /m b and a s {m h ). Corrections of order Aq CD /m 6 to T 12 have been calculated 
in PHTU] and applied to Qf s in [5]. 



2 rf 2 at next-to-leading order in QCD 



We decompose r 12 as 



12 



\ 2 T^CC 

A c 1 12 



^ 12 + ^u 1 12 



(3) 



with the CKM factors A, = V/ d Vib for 2 = u,c,t. The coefficients T^, a, b = u, c in Eq. (jHJ), 
which are computed from diagrams like those in Figure d are positive. We present the new 
NLO expressions for the coefficient in the appendix. has already been given at NLO 
in [7j, and r" 2 can be inferred by taking the limit z — > in V^ 2 . It is convenient to write 
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1 12 
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Here B = £>(/U 2 = m&) and .B^. = B s (fi 2 
the local AS = 2 operators Q and Q5: 



B +am ) \ t 
nib) parameterize the hadronic matrix elements of 



Q 

(B d \Q(u, 2 )\B d ) 
(B d \Q s (u, 2 )\B d ) 



?7m(! -l5)bq^(l -75)6, = 9(1 + 7s)6g(l + 7s)&, 



-/I ^1 



B(^), 



■X 



2 M 2 



The mass M Bd and decay constant /e d cancel from Eq. (j3J). I? and B' s depend on the scale \i 2 
and the renormalization scheme used in the computation of the matrix elements in Eq. JSJ). 
When combining values for B' s /B with our results for ci )2 , ai j2 and &i i2 below, one must verify 
that they correspond to the same scheme. Details on the renormalization scheme used by us 
can be found in Often the parameter B$ rather than B' s is chosen to parameterize (Qs)- 
As shown in Eq. (jJJJ), they differ by a factor involving MS masses. mb(m&) is smaller than the 
pole mass m& by roughly 0.4 GeV. 

For the evaluation of Eq. (j3J) we also need the SM prediction for Mi 2 : 

t2 \ 



B, 



/I.M 2 



Mr 



>B S (»2). (5) 



M 



12 



A? Gf 



* 12tt 2 

with the QCD factors ?7b = 0.55 [S] and 

6 S 0) = M/i)]" 6/23 



M BdVB B(^)b B ^ 2 )f 2 B M^S 



M. 



(6) 



1 + 



a s (/x) 5165 



b B (m b ) 



1.52 ±0.03. 



4tt 3174 

Note that results from lattice gauge theory are often quoted for the scale and scheme invariant 
parameter B = b B (fi2)B(fi 2 ) rather than B{m^) entering Eq. (j3J). 

We use the following input for the physical parameters (where m 8 = m^m;)): 

m b = (4.25 ± 0.08) GeV, m c = (1.30 ± 0.05) GeV, 
a s (M z ) = 0.118 ±0.003, m t = (167 ±5) GeV, 

B' s /B = 1.4 ±0.2, ml ow = (4.8 ± 0.2) GeV. (7) 
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The top mass mainly enters the result through Sim^/My/) in Eq. (jSJ, which evaluates to 
S(rrit/Myy) = 2.40 ± 0.11. In the power corrections a m , b m , c m the renormalization scheme is 
not fixed, because corrections of order a s jm^, are unknown. The expansion parameter of the 
HQE is the pole mass and we use ml° w = 4.8 ± 0.2 GeV (and m d = 7MeV) in a m , b m and c m . 
For the determination of 

a = a 1 + a 2 — + a m (8) 

and the analogously defined quantities b and c we take B' s /B = 1.4 ± 0.2, which covers the 
range of recent lattice computations [T2]. We estimate the accuracy of our calculation by 
computing the coefficients in two schemes for the quark masses (pole and MS), as explained 
in the appendix. Further we vary the renormalization scale /ii between one half and twice the 
b quark mass in the corresponding scheme. The result is shown in Figure |2] for the coefficient 
a, which is most relevant to af s : While the dependence on \x\ is small in both LO and NLO, 
the scheme dependence is huge in LO and reduced by roughly a factor of 4 in NLO. We quote 
our coefficients for the two schemes and add the errors from Eq. (JJJ), and the uncertainty from 
the /^-dependence in quadrature: 
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(9) 



In the case of a m , . . . , c m the difference between the LO and NLO columns stems solely from 
the QCD factor t)b- The reduction of the scheme dependence of a±, . . . , c 2 is evident from the 
comparison of the last two columns with the first two ones. 

Our final values for a, b, and c are at NLO (LO results in parentheses): 

a = 12.0 ±2.4 (14.7 ±6.7) 

b = 0.2 ±0.1 (0.6 ±0.5) 

c = -40.1 ±15.8 (-63.3 ±15.6) (10) 

They have been obtained by averaging the results in the pole scheme and the MS scheme 
for central values of the input parameters. The error from scheme dependence was taken to 



2 rf 2 at next-to-leading order in QCD 
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Figure 2: Dependence of a on the scale /ii. The solid (dashed) lines show the NLO (LO) 
results. 

be half the difference between the results in the two schemes. The errors quoted in Eq. ([10)1 
were obtained by combining in quadrature the latter error with the uncertainties in the MS 
scheme from scale dependence (fii), fn c , mb, fn t , ct s (Mz), B' s /B, and the 6-mass in the power 
corrections. 

In order to understand the size of the coefficients a, b, c at leading and next-to-leading 
order and the impact of various uncertainties, it is instructive to expand in the small parameter 
z = ml/ml ~ 0.1. The leading terms in this expansion behave as follows: 
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Here we have displayed the coefficients a i: bi and q separately, indicating the leading order 
terms and the NLO corrections. 

In the SM the CP asymmetry Of s does not depend on Cj, but only on a; and bi, on which we 
shall focus for the moment. Both a and b exhibit an interesting pattern of GIM suppression, 
which leads to a pronounced hierarchy among the different contributions. All of the coefficients 
of Ofg have to vanish as z — > 0. The dominant term is ai, while a 2 is suppressed by one, &i j2 even 
by two additional powers of z at LO. This strong hierarchy is alleviated at NLO, where the 
z 2 and z 3 terms receive corrections of order a s z and a s z 2 . Hence they are still parametrically 
smaller than a±, which remains the most important coefficient. As a consequence of this 
pattern, the coefficients &i )2 get larger relative corrections at NLO, but remain strongly 
suppressed in comparison to a\. This suppression is also not changed by the power corrections 
b m . Thus b has only a minor impact on af s . An additional welcome feature is the suppression 
of a 2 , which considerably reduces the dependence on the hadronic matrix elements B' s /B. 



LO[pole] 




NLO [pole] 

NLO[MS] 



LO[MS] 
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We emphasize that the dominant term a\ is free of hadronic uncertainties since the matrix 
element B in T 12 cancels against the identical quantity in M 12 . It can be seen from Eq. 
that power corrections to a are suppressed by an additional factor of z. As a result of all these 
properties, af s is quite accurately known in the SM, once the NLO QCD effects are taken into 
account. Note that the latter are important to eliminate the sizable scheme ambiguity of the 
leading order calculation. We remark that the a s zlnz term in a\ is peculiar to the choice of 
pole masses z = fn 2 ole /ml pole , which at one-loop order is equivalent to z = m^(m c ) /rrif(m&) . 
Expressing the results in terms of z — m 2 c {jn^) jrfl^ijn^) , the zhiz term is eliminated. As 
discussed in [TJ the absence of these terms holds to all orders in a s . Finally, at NLO the 
overall uncertainty in a and b comes predominantly from m c and from the residual scheme 
dependence. 

The situation is different for c, which is enhanced relative to a, b. Here sizable uncertainties 
are still present at NLO from the dependence on B' s /B, power corrections and, to a lesser 
extent, also from residual scale and scheme dependence. The parameter c enters the width 
difference Ar^ and, in general, the expression for cif s in the presence of new physics. In these 
cases one has larger theoretical uncertainties than in the SM analysis of <2f s . 



3 Phenomenology 

In the SM the CP asymmetry for the system reads 

Ti 2 
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12 
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(12) 



where a and b are given in Eq. (|10[). In terms of Wolfenstein parameters p and rj the CKM 
quantities in Eq. (|T2"|) are 
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Ri 



m 



(13) 



(14) 



where P = arg(— A t /A c ) and R t = y(l — p) 2 + ff are one angle and one side of the usual 
unitarity triangle. 

A future measurement of af s will allow us to constrain p and f] within the SM using the 
theoretical values for a and b. This is illustrated in Figure El 
Using Eq. (HUJ) and [HI 



we predict for af s in the SM 



i?t = 0.91 ±0.05, p= (22.4 ± 1.4)° 

af s = -(5.0 ± 1.1) x 10~ 4 



(15) 
(16) 
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Figure 3: Constraints in the (p, fj) plane implied by given values of the CP asymmetry af s . 
The area between the solid pair of curves on the right represents the theoretical uncertainty 
at NLO, assuming af s = — 1CT 3 . Similarly, the curves on the left indicate the uncertainty for 
af s = — 5 x 10" 4 both at NLO (solid) and at LO (dashed). The currently favoured solution for 
the unitarity triangle is also shown. 



This result is entirely dominated by the a-term in Eq. (|12j) since the small contribution from 
b is further suppressed by its CKM coefficient, which is small for standard CKM parameters. 

Our results can also be applied to the case of B s mesons, where Eq. (jl2j) holds with 
obvious replacements. Here the term proportional to b is strongly CKM suppressed and can 
be neglected. SU(3) breaking in a is negligible as well and the result in Eq. (fTUj) may be used. 
We then find (V us = 0.222) 

a s k = a\V us \ 2 R t smp x 10~ 4 = (0.21 ± 0.04) x 10" 4 (17) 

The width difference in the Bd system is given by ATd/AMd = — Re(ri2/Mi2). The real 
part of T12/M12 can be found using Eqs. (jU), (fTUj) . (fT3*j) and (JT3J). It turns out that for the 
parameters in Eq. (fl3j) the c-term yields the full result to within about 2%. In view of the 
large uncertainty of c, the contributions from a and b can be safely neglected. We then obtain 
the SM prediction 

Ar d , ln , , 3 Ar d 



AM d 



(4.0 ± 1.6) x 10"% = (3.0 ± 1.2) x 10" 



where the second expression follows with the experimental value AM^/r^ = 0.755. This result 
for ATd/Td is in agreement with [T|ll()j. To the extent that SU(3) breaking in the ratio of bag 
factors B' s /B can be neglected, the number for AT / AM in Eq. ()18j) applies to the B s system 
as well. 

The effects of new physics in M i2 on af s have been discussed in [5] . If magnitude and phase 
of M12 are parameterized as 

Mi 2 = r 2 d e 2W «Mf 2 M (19) 
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one obtains jH] 



\M 12 J SM rj \M 12 J SM r\ V ; 



Since the real part of r 12 /M 12 in the SM is much larger than the imaginary part, af s is 
particularly sensitive to new physics. In this more general context our results can also be 
used. However, it has to be kept in mind that the SM analysis leading to Eq. (|TK|) may no 
longer be true in the presence of new physics and the determination of CKM quantities then 
needs to be modified. 



To summarize, we have computed the CP violating observables aj s at next-to-leading order 
in QCD. We include the effect of penguin operators in the weak Hamiltonian and the power 
corrections of relative order kQ CD jm h . Our SM predictions are given in Eqs. (fTBj) and (JT7J). 
We emphasize that within the heavy-quark expansion the a| can be reliably computed in the 
SM as functions of CKM parameters. A crucial element is the small sensitivity to hadronic 
parameters, which enter only as the ratio B' s /B and only with a suppression factor of z = 
(me/mi,) 2 . After including the NLO corrections, the theoretical error on a| s is reduced to about 
20%. This is largely due to a reduction of the scheme ambiguity in the definition of quark 
masses by a factor of 4 in comparison with the LO result. The remaining uncertainty is larger 
for AT q . The result at next-to-leading order in QCD is given in Eq. (fTHj) . The measurement 
of a]?, is possible using suitable flavour-specific decay modes of neutral B mesons. If it can be 
performed with sufficient accuracy, it will provide a significant test of the Standard Model. 
The large sensitivity of a q {s to new physics is reinforced by the improved theoretical analysis 
presented here. 



Note added 

The topic of this paper has also been addressed by Ciuchini et al. [H], who pointed out an 
error in an earlier preprint version of this paper. Our analytical results in Eq. (|25|) now agree 
with those in Eqs. (43-45) of We thank the authors of ^3] for clarifying communication. 
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A NLO coefficients 

Here we collect more detailed results for the coefficients in Eq. Q. The HQE expresses T 
for the Bd system as 



ah 
12 



Tttfc 

12 



fl M Bd 
24vr JBd d 



F ab (z) + P ab (z 



-B 

3 



Ff(z) + P°, 



+ r 



ab 

12,l/m 6 



.(21) 



A NLO coefficients 



9 



The short-distance coefficients F (z) contain the contributions from the AB = 1 current- 
current operators Q\ and Q 2 - The NLO results for F cc (z) and Fg c (z) have been derived 
in P, where these coefficients are called F(z) and F s (z), respectively. Further F uu = F cc (0) 
and Fg u 



Fg c (0). The coefficients P(z) and Ps(z) contain the contributions from penguin 



operators. They come with small coefficients, which simplifies the NLO calculation 
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with x, 



fi/rrih and an analogous notation for F^. The AB = n operators, n = 1,2, are 
defined at the scale /x n = Oim^). The dependence of F^j c on /ix diminishes order-by-order in 
a s . 

Throughout this paper we use the same operator definitions and renormalization schemes 
as in [Zj, with one important addition: In af s the renormalization scheme of the quark masses 
is an important issue and we choose two different schemes for the computation of the a^, bi, 
Ci in Eq. (@J). For both schemes we take the MS masses m c (m c ) and m^rrTf,) as the basic 
input. In the first scheme (pole scheme) we express the observables in terms of nib = m^poie = 
m&(l + 4a s (m;,)/37r), using the one-loop relation between pole- and MS-quark mass. In this 
scheme we define the variable z as z = (m c (m c ) /m^mf,)) 2 , which to one-loop order is equivalent 
to the ratio of pole masses squared. In the second scheme (MS scheme) we take m& = mb{jnb) 
and replace z by ~z = (m c (rrib) /mb(m&)) 2 , where both running masses are defined at the 



scale mft. The results below for the functions F { 
corresponding functions F 
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are valid in the pole scheme. The 



in the MS scheme are obtained via the relation 
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The coefficients read: 
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Li 2 (z) + 



ln(l - z) ln(z) 



+ 



4(1- z) 2 (1 + z) (-l + 13z + 3z 2 
3~2 

-16(l-z) 2 (l + 2z)] ln(x Ml ) + 
8 (1-z) (1 + 2^) 



ln(l-z) + 



4z (2-3z + 18z 2 -3z 3 



32 (1-z) 2 (l + 2z) 



ln (^ 2 ) + 



2 28 (1 - 2) (-5 - 82 + 192 2 ) 

7T + 



ln(z) + 



(25) 



In terms of the function P(z) used in j?j the penguin coefficients in Eq. ()21|) read P cc (z) 
P{z), P uu = P(0) and 



P| c (z) 



P s (z) + P s (0) 



8AP 1 



(26) 



with 

AP UC 



a s {m) 1 - (1 + 2*^1 - 4z 



47T 



C 2 ( Ml ) 



In ,2 - (1 + 2z)Vl -4zlna - Az] (27) 



and cr = (1 — \/l — Az) / (1 + a/1 — 4z). AP MC is of order and numerically negligible. 

The power corrections E^ y were first obtained for ab = cc, uu in and for a& = uc 
in [TH]. We have re-computed the case ab = uc here, confirming the results of [10]. in the 
notation of we find ((. . .} = (B\ . . . \B)) 



1 12,l/m 6 



2AnM B 

1 + z + z 2 



(1-z) 5 



(1 + 2z)K 2 (R ) - 2(1 + 2z)(K 1 {R 1 ) + K 2 (Pi)) 



-2 



1 - z 



(K^Rz) + K 2 (R 2 )) 



12z 2 
1 - z 



(K^Rs) + K 2 (R 3 )) 



(21 
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